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Introduction 

Adaptive-grid methods are an important subject of study in computational fluid dynamics because of their 
potential for improving the efficiency and accuracy of numerical methods. Various papers on grid adaptation 
have been presented in recent years, but few methods have been applied in practical applications. 

In this paper, a practical solution-adaptive-grid method utilizing a tension and torsion spring analogy is 
proposed for multidimensional fluid-flow problems. The tension spring, which connects adjacent grid points 
to each other, controls grid spacings so that clustering is realized in regions containing shock waves and 
shear layers. The torsion spring, which is attached to each grid node, controls inclinations of coordinate lines 
and prevents excessive grid skewness. A marching procedure is used that results in a simple tridiagonal 
system of equations at each coordinate line to determine grid-point distribution. Multidirectional grid 
adaptation to the flow solution is achieved by successive applications of one-directional adaptation. Examples 
of applications for axisymmetric afterbody flow fields and two-dimensional transonic airfoil flow fields are 
shown. 

Grid Adaptation in One Coordinate Direction 

For simplicity of illustration, consider grid adaptation to a flow field in which grid points are free to move 
along each q-coordinate line whose configuration is fixed. Let a grid point A in figure 1 be connected to its 
adjacent points, B and C, by tension springs whose spring constants are and K tJ . To distribute grid 

points along the q< -coordinate line in proportion to the gradient of selected flow properties, the relationship 
between the spring constant K and the gradient of the dependent variable / 

K %J = 1 + — *»,j) (1) 

is used, where C\ is a constant and s,j is the arc length calculated from point (i,l) along the q, -coordinate. 

Using equation (1), the distribution of grid points along the q, -coordinate line, namely, new values of s tiJ , 
is determined by 

' ®».j) — ~ = 0 (2) 

The idea of a spring analogy represented in equation (1) was introduced by Gnoffo(ref. 1). In his model, 
the distribution of points along the q, -coordinate is determined only by the gradients of flow properties 
along that coordinate line and is not affected by the distribution on adjacent q-coordinates, q,-i and q,+i. 
This can lead to excessive skewness of grid lines, especially when applied to complex flow fields, and this lack 
of control of grid inclination makes it difficult to extend the scheme to more than one family of coordinates. 

A force to control inclinations of ^-coordinates in addition to that of grid spacings on q-coordinates will 
correct this deficiency. This control force can be given by consideri ng to rsion springs attached to nodes 
along the q,_ x line. The torsion spring enforces the inclination of line DA to that of a reference line. If the 
spring constant of the torsion spring is denoted by IT, a mathematical statement of the force is 


Ftortion — [^DA ~ <f>) (3) 

where 6da is a inclinatio n of l ine DA and <f> the inclination of the reference line. The reference line can be 
chosen as an extension of FD to avoid kinks in the (-line at point D, as a line normal to the q,-coordinate 
to make the grid quasi-orthogonal, or as a streamline, and so forth. In practical calculations, a combination 
of these reference lines is used. The torsion spring constant H can be prescribed for each coordinate line. 
A balance equation for the complete spring system is 

— 8 t,j) ~ K*,)— l(*M — — -^»-l J ~ fit-l,]) = 0 


(4) 



7b facilitate solutions to equation (4), the third term is rewritten 

~ => (5) 

where is arc length to the intersection of reference line DA' with t he q, -coordinate as depicted in figure 
2. The term is set equal to 3 t -ij divided by length of DA 1 . Finally, equation (4) reduces to the 

following equation: 

■ffiij— l 8 *,)— 1 ” i^iil + Ktj—l + (6) 

This is a tridiagonal system of equations for 8 tt) and can be readily solved. 

In this analysis only the torsion force on the upstream side (q,_i) influences the distribution at q,. This 
permits simple marching schemes to be used, without any loss of generality, and contributes to the simplicity 
and robustness of the method. Conversely, if the influence from both sides (q,_i and tfr+i) is considered 
simultaneously, the computational effort is increased considerably without any additional benefit. Note, 
too, that the downstream influence, (q, + i), could be used instead of the upstream, without any additional 
complexity or loss of generality. 

Grid Adaptation in Multicoordinate Direction 

A model using tension and torsion springs for a two-directional adaptation can be depicted as in figure 
3(a). Since each grid point is connected to its four adjacent points, the procedure for the grid movements of 
this model is more complicated and requires more computational effort. To minimize this complexity, a split 
model, which is a combination of one-directional adaptation, is used(fig. 3(b)). Grid movement is achieved 
by successive applications of the one-directional adaptation method. This is analogous to ADI schemes for 
partial differential equations. It is not necessary to achieve convergence in this adaptation procedure and, 
in fact, one iteration is sufficient in practical applications of the scheme. 

The solution field is interpolated onto the newly adapted grid, using second-order, one-dimensional La- 
grange interpolation after each one-directional adaptation. Before this interpolation, it is possible to add 
or delete grid points at the users discretion, thus enhancing the method without any loss of accuracy or in- 
creased complexity. The successive application of one-directional adaptation also enhances the applicability 
of the method, in that unidirectional adaptation can be used in regions where gradients are large in only 
one direction. Extension of this scheme to three-dimensions is straightforward. 

Results 

Shown in figures 4-8 are examples of applications of the method to axisymmetric, plume flow fields. 
These complex flow fields exhibit oblique shocks, barrel shocks, and slip surfaces, where grid points should 
be clustered for adequate resolution. The locations of these discontinuities are not known apriori, and a 
solution-adaptive grid is particularly useful in realizing an efficient and accurate simulation. 

The initial grid(fig. 4) was generated by an algebraic method and the flow field determined using a 
code(ref. 2) for the thin-layer Navier-Stokes equations. The free-stream Mach number is 2.01, the jet-exit 
Mach number is 2.5, and the static pressure ratio of the exhaust jet to the free-stream is 1. The density 
gradient was chosen as a reference variable for the grid clustering, and figure 5 shows a solution-adapted 
grid which clearly has clustered points to the oblique shock, the slip surface, and the barrel shock. Shown 
in figure 6 are computed density contours using the grid shown in figure 5. 

Shown in figure 7 is an adapted grid for a high-jet-pressure case (pressure ratio of 6), and density contours 
are shown in figure 8. The same initial grid (fig. 4) used in the previous case was used here. 

The grids were adapted to the flow-field solutions periodically (typically two or three times) during the 
course of reaching a steady-state solution. In this way, the time required for adaptation is a negligible 
fraction of the total time required for flow-field solution. 

Examples of the application to transonic flow fields past a NACA0012 airfoil are shown in figures 9-12. 
The free-stream Mach number is 0.8 and the angle of attack is 1.25°. Two-dimensional Eulef equations 
were solved, using the code ARC2D(refe. 3,4). Figure 9 shows an initial O-grid generated by an algebraic 
method. Grid-point distributions along ^-coordinate lines, which are parallel to the airfoil surface, were 
adapted to the density gradients in the flow-field solution. Shown in figure 10 is an adapted grid showing 
appropriate clustering of grid points to the shocks. Computed Mach contours obtained using the adapted 
grid are shown in figure 11. Figure 12 shows a comparison of pressure coefficients obtained using nonadapted 
and adapted grids. This figure shows that shock waves are crisply resolved using the adapted grid. 

An example for supersonic flow is shown in figures 13-15. The free-stream Mach number is 1.2 and the 
angle of attack is 1.25°. In this case, both £ and q -coordinate lines were adapted to the density gradients. 



Summary 

The principal features of the proposed adaptive-grid method are as follows: 

1. It is a simple concept consisting of a tension and torsion spring analogy. The combination of these springs 
produces a suitable adaptive grid without excessive grid skewness. 

2. A marching-type of calculation procedure minimizes the computation time. 

3. The split-solution procedure for multidirectional adaptation is simple and practical. 

4. The method can be applied independently to selected parts of the entire grid. 

5. The controllability of grid inclinations with the torsion spring makes it possible to generate a nearly 
orthogonal adaptive grid. Also, the inclinations of grid lines near boundary can be specified arbitrarily. 

These features make the proposed adaptive-grid method practical and robust, and enhance its applica- 
bility. The method can be applied to two- and three-dimensional fluid-flow problems without any difficulty 
and requires little computational time and effort. 
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Fig.l Schematic of adaptive-grid algorithm 
with tension and torsion spring analogy. 





Fig.2 Notations for torsion spring analogy. 
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Split model. 


Fig.3 Multidirectional adaptation. 
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(a) Complete configuration. 


(b) Base region detail. 


Fig.4 Initial grid for afterbody flow field. 


!■ 

111 1 1 

III !■ 
iiniiilll 

. M int" iiillli 

!===:: !!!!S!!!l!!!| 


mill 

■nil 

Imii 

lanii 


lESS;; ! -j;:!!'! 


llllllll 
I III III 

IBM 








Fig.5 Adapted grid for afterbody flow field: 
Moo = 2.01, My = 2.5 : Pj/Poo = 1.0. 


Fig.6 Computed density contours with adapted grid: 
Moo — 2.01, Mj — 2.5 . Pj /Poo — 1.0. 
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Fig.7 Adapted grid for afterbody flow field: 
Moo = 2.01, My =25: Pj/P^ = 6.0. 


Fig.8 Computed density contours with adapted Q rid: 
Moo = 2.01, My = 2.5 : Pj/Poo = 6.0. 
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